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Abstract 

The statistics of the curvature quanta generated during a stage of inflationary expan- 
sion is used to derive a count response model for the large-scale phonons determining, in 
the concordance lore, the warmer and the cooler spots of the large-scale temperature in- 
homogeneities. The multiplicity distributions for the counting statistics are shown to be 
generically overdispersed in comparison with conventional Poissonian regressions. The gen- 
eralized count response model deduced hereunder accommodates an excess of correlations 
in the regime of high multiplicities and prompts dedicated analyses with forthcoming data 
collected by instruments of high angular resolution and high sensitivity to temperature vari- 
ations per pixel. 
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According to the conventional lore of structure formation, Cosmic Microwave Background 
(CMB) observations provide an image of quantum fluctuations blown up to the size of 
the Universe. Curvature perturbations are believed to originate as quantum fluctuations 
at least in the framework of the standard ACDM paradig where scalar modes of the 
geometry are the sole source of temperature inhomogeneities because of the strict absence 
of tensor modes. In this simplified scenario (yet consistent with the three independent 
cosmological data sets [H [2j [3]) the curvature perturbations are quantized in terms of a 
collection of scalar phonons and are themselves proportional, via the Sachs- Wolfe effect, to 
the temperature inhomogeneities. The latter statement holds for large angular scales (i.e. 
•d > 6 deg) corresponding, in practice, to the region of the so-called Sachs- Wolfe plateau, i.e. 
multipoles i < 30. This occurrence is often dubbed by saying that CMB maps at the largest 
scales are a faithful impression of quantum fluctuations as mentioned in the first sentence 
of this paragraph. For smaller angular scales (i.e., approximately, 30arcmin < d < 6 deg) 
the curvature perturbations are still the source of temperature anisotropies which enter the 
regime of acoustic oscillations. For even larger multipoles (i.e. i > is with is — 920) diffusive 
(Silk) damping dominates. 

The effective action obeyed by the primordial phonons during a conventional stage of 
inflationary expansion can be written, in terms of 7Z, i.e. the curvature perturbations on 
comoving orthogonal hypersurfaces which have a gauge-invariant meaning without being 
necessarily connected to curvature perturbations in all the coordinate systems different from 
the comoving orthogonal one (see, for instance, [U [5] and references therein): 



where <fi denotes the (single) inflaton, a is the scale factor and 7i = (In a)'; the prime denotes 
a derivation with respect to the conformal time coordinate r. In Eq. ([1]) the geometry is 
assumed to be spatially flat, as suggested by the position of the first acoustic oscillation 
and by CMB data as a whole (see, e. g. [I]). In the minimal ACDM scenario (when the 
non-adiabatic pressure fluctuations are strictly vanishing) the curvature perturbations are 
approximately constant, i.e. 7Z' ~ for wavelengths shorter than Hubble radius at each 
corresponding time (see, for instance, [6]). In particular, at photon decoupling, the intensity 
fluctuations of the radiation field (i.e. the warmer and the cooler regions in the CMB sky) 
are given, in real space, as 



under the approximation of sudden decoupling (i.e. assuming that the visibility function 
is a narrow Gaussian centered at Td ec )- Following the tenets of canonical quantization, 

2 In the concordance model (often dubbed ACDM where A stands for the dark energy component and 
CDM for the dark matter component) single field inflationary models and standard thermal history are 
always assumed implicitly. 
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the field variables can be promoted to quantum mechanical operators obeying equal-time 
commutation relations, i.e. 



K(x, t) 



(2tt) 3 / 2 VV 



-ik-x 



z(t)\[W 
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where d^ and a) p are the annihilation and creation operators of curvature inhomogeneities 
and at] = ^ • per each Fourier mode the Hamiltonian operator describing the quantized 
curvature perturbations can be written as: 



H % =2k lC {k) + 2 



A*(r)/C_(fc) + A(r)/C+(fc) 



(4) 



where 2 A = iz'/z; the operators /C±(/c) and JCo(k) obey the commutation relations of the 
SU (1,1) Lie algebra [ZHSJ: 



From Eq. (j4]) the multiparticle final state is given by 



at dg + d_gd^ 



1%) = s(^) s(c fc )|o ff o_^, |*> = n 1%), 

k 

where the two unitary operators ^((fi k ) and S(Cfe) are 

S(<p fc ) = exp [-2«<p fc /C (fc)], S(Cfc) = exp [QJC_(k) - ( h }C + (k)], 



(5) 
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with ( k = r k e %lk and a k = (2ip k — j k ); the time evolution of the variables r k (r), <fk( T ) and 
afc(r) is given by 



= 2iA cos a&, ip' k = k — 2i\ tanh sin a^, 



a' k = 2k- 4iX Sin i Z k ■ (8) 
k tanh2r fc v ' 



In Eq. ([6]) the initial state is the vacuum. While the latter assumption can be relaxed, it is 
usually invoked by tacitly assuming that the total number of inflationary efolds exceeds the 
minimal amount required to solve the standard problems of big-bang cosmology (see, e.g. 
[I], first reference). 

The customary exercise would now be to compute the power spectrum in terms of the 
parameters of the underlying inflationary model encoded in the function z(t) defined in 
Eq. ([Tj). A qualitatively different class of questions concerns instead the determination 
of the multiplicity distribution of the curvature quanta. Owing to the 577(1,1) group 
structure of Eq. (j^D, the squeezing operator of £(Cfc) defined in Eq. (jUJ) can be fac- 
torized as £(( fe ) = A+(Ck) ^4o(Cfe)^-(Cfe) where A(Cfc) = exp [-2 In (coshr fc )/C (fc)] and 
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>A±(Ck) = exp [=pe ±i7fe tanhr^ /C-|-(fc)]. From Eq. ((H]) the form of the density matrix rele- 
vant for the forthcoming discussions is 

1 oo oo 

ft = E E e-^^^CtanhrO^lnj n_ % ){m_ % (9) 

cusn ffe ng=0mg=0 

whose diagonal elements define the multiplicity distribution: 

— n k 

n/ 

*W = It Pn V Pn % {n % ) = (1+ ^ rs+ l . (10) 

accounting for the way curvature quanta of each Fourier mode (i.e. nr) are distributed as 
a function of their mean value per each Fourier mode (i.e. = sinh 2 r^,). Equation (flQ|) is 
a Bose-Einstein distribution but the average number of curvature quanta has no relation 
with the standard Bose-Einstein occupation number. The same situation occurs usually in 
quantum optics for chaotic (i.e. white) light |9] HD] where photons distributed as in Eq. ( TlQjl 
for each mode of the radiation field can be produced by sources in which atoms are kept at an 
excitation level higher than that in thermal equilibrium. The way the off-diagonal elements 
of Eq. (jUJ) behave is dictated by the phases a^. While it is not strictly essential to get rid of 
the off-diagonal elements for the forthcoming arguments, it is nonetheless plausible that, by 
averaging over Eq. the density matrix can be reduced (i.e. p*? d = ^ J 27r dotkf)^) by 
only keeping the diagonal terms. The density operator can be written, with the conventional 
shorthand notation, 

P = E P{n,} |{n 5 }><to}|, |{nj) = In^K-JK-;)... (11) 

where, P{ n -} is given by Eq. ( ITOi) and the ellipses stand for the direct product of all the 
modes of the field. According to Eq. (J2j) (and in the hypothesis that curvature quanta are 
the sole source of temperature inhomogeneities) , the distribution P(n) of the total number of 
phonons n = Ylt n k m ust reflect the distribution of the warmer and cooler regions. If there are 
supplementary sources of temperature inhomogeneities (e. g. tensor modes of the geometry) 
the present discussion can be appropriately modified. The multiplicity distribution P{n) 
accounts for the way the total number of phonons n is distributed as a function of its mean 
value; P(n) can be very different from P{ n ^}- Denoting with p{{n}) the joint probability 
distribution of the set of phonon occupation numbers {n} of the field, we shall have that 

p(W1 1 (i + %)(i' + V%)v (12) 

For any mode for which = 0, the corresponding factor must be interpreted as 5n £ o- I n 
the following we shall suppose, quite generally, that only a subset consisting of e modes of 
the field is actually occupied and we shall restrict the attention to this subset of modes. 
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If n is the total number of phonons and P(n) is the multiplicity distribution of n, then 
P{m) = E{n}P(W)^mn' In quantum optics an analog of the multiplicity distribution P(m) 
describes the statistical properties of (unpolarized) chaotic light beams [9]. The evaluation 
of P(m) can be in general difficult but it becomes easy in the physical case when the average 
occupation number of all the e occupied modes become equao; in this case from Eq. ( TTTT) 
and (TT2"j) the joint probability distribution of the occupied modes becomes: 

p({n}) = — r^— , n = y^nr = eni:. (13) 

IJ (1 + n/e) e (l + e/n) n ' £ 

Every non- vanishing term in the summation Pirn) = J2{ n } p({ n })$mn has the same value 
and the required probability P(m) is simply p{{n}) given by Eq. ([TBI multiplied by a well 
known combinatorial factor accounting for the way the n phonons are distributed among the 
e modes: 

T(n + e) f n \ n { e Y , , 



(15) 



r(e)T(n + 1) \n + ej Vn + e, 
The cumulant generating function [TT] associated with Eq. ( |T4"1) is given by 



C(s, n, e) = In 



n,e 



— e In 



! + (!-«)- 
e 



Equations (U4j) and ( 1151) define the count response model implied by the physical nature of 
the source which are the (quantized) curvature perturbations. From the cumulant generating 
function all the cumulant moments can be obtained to an arbitrary order but they are all 
function of n and of e. The variance D 2 = (n 2 ) — (n) 2 is then given by D 2 = n + n 2 /e. In 
the limit e — > 1 the Bose-Einstein distribution is recovered; for generic e the count response 
model defined by Eqs. (fT4"l) and ( TT5l) falls into the class of negative binomial regressions 
which arise, in rather general terms, in all those discrete counts where correlations lead to 
asymmetric distributions with a degree of correlation larger than in conventional Poissonian 
counting (see [12] for an introduction to statistical models of count response data). 

The count response model of Eqs. ( JT41 and ( 1T51) has been deduced in a top-down approach 
by looking at the statistical properties of the counting distribution of primordial phonons as 
they arise in the minimal ACDM model. A complementary avenue will now be taken with 
the purpose of deriving the concept of multiplicity distribution in a bottom-up perspective. 
To study the multiplicity distribution of the spots we can, for instance, fix a threshold in 
the brightness perturbations such as 

(nun), _ na QaK f -An \ 1/2 { T, /0 



|A m j |= 26.865 -^-^) UK, (16) 

1 \2.43 x 10~ 9 / V2.725K/ P 



3 In the case of a thermal light beam which is either fully polarized or fully unpolarized the use of a 
rectangular spectral density is an excellent approximation in the derivation of the photocounting statistics 



which is also experimentally accessible [9]. 
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where An is the amplitude of the power spectrum of curvature perturbations following 
from the WMAP 7 data at the conventional pivot scale k p = 0.002 Mpc -1 . The value of 
Eq. (fl6|) comes from the large-scale ($ > 6deg) plateau; different and more refined ways 
of fixing the threshold can be suggested but this aspect is immaterial for the forthcoming 
considerations @. Given the threshold (FTO]) we can therefore ask (or predict) how many hot 
(or cold) spots are present in different angular intervals starting from a $ m m (connected, 
for instance, with the resolution of the instrument) up to a i9 max - The resolution^ will also 
determine ultimately the nature of the partition and the number of classes of the histogram. 
The same procedure will then to be used for different temperature thresholds. It is plausible 
to assume, in a bottom-up approach, that the distribution of spots in excess with respect 
to a given (progressively increasing) threshold is given by a real Gaussian variable /($). In 
a Poissonian count response model the probability that one spot is found between i? and 
i? + d$ will be P(l, d-d) = d$ where /x is a constant and [11] 

QW = t/ f 2 (#)d#. (17) 



5 

The probability of finding n spots in the interval d + At9] is a Poisson distribution in n 



F(„,0,0) = I 



r-d+A-d 1 n 

H / Q(d)d§ exp -u / Q(d') d& 

. Jd J L J$ 



(18) 



It would be debatable to conclude that the multiplicity distribution is Poissonian since what 
is potentially measurable is not P(n, Ai?) but rather p n (A$) = (P(n, Ai?)), i.e. the 
distribution P(n, A$) expressed as a function of x = // +Atf Q(^) dt^ and averaged over 
the ensemble of x. The mean and the variance of x can be physically estimated as (x) = Q At9 

2 

and ) - (x) 2 = Q (At?) t9 c where $ c is a typical scale possibly related to the angular 

resolution to some other coarsening angle. The question is therefore if it exists a probability 
distribution of x (containing at least two parameters) and reproducing the results already 
obtained within our top-down approach. The answer to the latter question is provided by 
the Gamma distribution ITT] 



e~ Xx 1 A$ 

Pr(A, v\ x) dx = (XxY^—^diXx), A = =— , v = — . (19) 



The distribution of Eq. ( ITHl) must then be averaged over the ensemble of x 

1 r°° ( 1 A 



Vn (M) = - Pv (=-,—;x)e-^^x) n dx. (20) 



4 We are here focussing on hot spots but the same discussion can be conducted for cooler regions by 
considering the distribution of spots below a given (progressively decreasing) threshold. This the meaning 
of the absolute value in Eq. (Til)]) . 

The WMAP experiment in his five channels has a resolution which varies between about 1 deg (for the 
low frequency channels) to less than 0.2 deg for the high-frequency channels. Th Planck explorer experiment, 
in his nine frequency channels has a resolution which varies between few arcminutes (in the high frequency 
instrument) and about half a degree of the low frequency channels. 
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The explicit result of the integral indicated in Eq. (12"U|) is exactly given by Eq. (|14p with 
n = fiQAv) and e = (Ai?)/i? c . The Poissonian counting is recovered from Eq. (fl5|) : in the 
limit e — > oo the variance tends to the mean value (i.e. D 2 — > n) and the cumulant generating 
function of Eq. (fT5|) tends to the Poissonian limit, i.e. C(s,n, e) — >■ n(s — 1). These two 
occurrences are sufficient to infer that, for e — > oo, P n (n,e) —> n n exp [— n]/n!. A posteriori 
we can say that the negative binomial counting deduced in Eq. (1141) (and, indirectly, in Eq. 
( f20|) ) is, at once, more general and more physical than the Poissonian counting. 

The probability of finding n spots in a given angular interval must depend on the thresh- 
old. By varying the threshold of Eq. ( Fl6l) as |Aj mm ' l | — > |Aj |?7 (where rj > 1) the values of 
n and e will depend upon 77 (in particular we can intuitively expect that n will decrease with 
rj). The probability generating function V of the negative binomial distribution satisfies 

dV 

~i; = - g{v > v) > (21) 

SP, V) = --^V\nV + ~ (-)V[1 - V x '% (22) 
ear] n dr]\e J 

which has the form of a reverse Kolmogorov equation and where V(n, e) = J2™=o Pn(n, e)- It 
is tempting to interpret rj as a continuous evolution parameter of an appropriate branching 
process [H]. In this case we should probably impose the boundary conditions ^(1, 77) = 
and G(V, 77) = (for V — > 0) and ask that Q(P, rj) factorizes as the product of a function of 
V and of a function of rj: in this way the branching process will be stationary in 77 and Eqs. 
(T2T1) ( T22T) imply that one (or both) the conditions are satisfied 

1 de e 2 d fn\ „ , „ m . , . 

-^r Cl nT,{lh ™V = c 2 V(l-V>") (23) 

where C\ and C2 are two numerical constants. This means that, assuming that the stochastic 
processes is stationary the 77-dependence of n determines also the //-dependence of e. For 
instance, from the first relation of Eq. (I23p we have that l/e(r]) + lne(r?) = ai + h\ In 72(77). 

The modest but novel purpose of the present paper could be summarized by saying that 
to count spots in a physically meaningful way we have to understand which is the appro- 
priate count response model. It has then been suggested how the multiplicity distribution 
of the total number of phonons is connected to the multiplicity distribution of the CMB 
spots in different angular intervals and for different temperature thresholds. These physical 
considerations pin down a specific count response model which is overdispersed in compari- 
son with a naive Poissonian counting. While the present considerations can be generalized 
to related frameworks it seems interesting to pursue dedicated analyses aimed at measuring 
the multiplicity distribution of CMB spots. The experimental scrutiny will therefore have 
to infer the likely values of e and n for different thresholds rj and for different angular inter- 
vals (and, presumably, also for different frequency channels of the instrument). The latter 
program assumes that the resulting CMB signal used to extract the multiplicity distribution 
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will be free of the contamination of foreground (e.g. point-like) sources, a potentially severe 
problem whose discussion and is beyond the scopes of the present paper. 
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